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INTRODUCTION 

This  report  contains  results  of  numerical  tests  on  a method  of  moments 
algorithm  for  the  solution  of  thin  wire  radiation  and  scattering  problems. 
The  algorithm  has  convergence  behavior  similar  to  a highly  overdetermined 
field  point  matching  scheme.  It  has  been  assumed,  therefore,  that  informa- 
tion about  the  range  of  applicability  of  overdetermined  wire  formulations 
can  be  inferred  from  the  results  of  these  tests.  Unfortunately,  in  the 
tests  performed  so  far,  the  algorithm  has  proven  to  be  generally  inferior 

i 

I 

to  the  other  formulations  tested. 

The  report  is  divided  into  three  sections.  The  algorithm,  called  MBCRE 
(for  Minimum  Boundary  Condition  Residual  Error)  is  described  in  Section  I, 
Section  11  explains  the  test  result  tables  and  also  the  connection  between 
the  MBCRE  algorithm  and  overdetermined  point  matching  schemes.  Section  III 

I is  a listing  of  the  computer  program  used  for  the  tests. 

« 

I 
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Section  I. 

The  integral  equation  relating  the  current  distribution  on  a thin  wire 
to  the  tangential  component  of  some  arbitrary  impressed  electromagnetic  fieid 


where  z and  z'  are  coordinates  of  distance  along  the  wire;  E^(z)  is  the  known 

tangential  component  of  the  Impressed  electric  field  along  the  cylindrical 

wire  surface;  I(z')  is  the  (unknown)  current  along  the  wire  axis;  v*  is  t tie 

wire  length;  e,U  are  the  material  parameters  (electric  permittivity  and 

magnetic  permeability)  of  the  medium  in  which  the  wire  is  imbedded;  . is  the 

radian  frequency;  a is  the  wire  radius;  and  k (=  u kp^)  is  the  wave  number. 

The  algorithm  discussed  in  this  report  for  the  solution  of  (1)  assumes 

a common  moments  approximation  1 (z')  to  the  axial  current  Kz")  of  the  form 

a 

N 

I (zO  = I I.f.(z')  (2) 

a 1=1  ^ ^ 

where  the  f^(z')  are  piecewise  sinusoidal  expansion  functions  given  by 


2 


.1 


f . ( Z*) 
1 


sin  k(z'  - z 
sin  k(z^  - z^) 

sin  k(z^  - z*) 
sin  k(z^  - z^) 


Zj^  < z < z^ 


Z2  < z < z,^ 


(3) 


elsewhere 


and  the  are  a set  of  unknown  amplitudes.  (See  Fig.  1) 

When  the  approximation  I (z)  is  inserted  in  (1)  it  will  produce  a tan- 

H 

gential  field  along  the  wire  differing  from  tlie  forcing  function  E (z)  by  a 

z 

residual  error  r(z), 

p 

.2 


r(z)  = E (z)  G(z,z')I  (z')dz' 

z COG  a 


dz 


(4) 


It  is  desired  to  find  the  set  of  for  which  the  integrated  mean  square 


error  is  minimum,  i.e. 


£ 


r*(z)r(z)dz  (*)  denotes  complex  conjugate 


(5) 


From  equations  (3),  (4)  and  (5)  it  can  be  seen  that  e.^.  is  a positive  definite 
real  quadratic  function  of  the  with  a unique  minimum  at  the  point  where: 


9x, 


= 0 


!!i 


= 0 


(6) 


Applying  condition  (6)  to  equation  (5)  gives: 


8x 


(7a) 


3 


1 (r*(E) 

"''k 


+ 


ar*(z) 


r(z))ciz  = 0 


k = 1,2 N. 


rhe  derivatives  in  (7)  have  the  form 


9r  (z) 
3r  (z) 


-hk(z) 


a -ih,  (z) 


;ihf(z) 

DXk  k 


i,k(.)=  ^ I (;(z,z')f^(z')dz' 

dz  ■' 


0 


(7b) 


(8) 


Inserting  these  back  in  (7)  and  simplifying  the  result  gives 

-2  Re[  i r(z)hj^(z)dzl  = 0 
' 0 


2i 


or 


0 


? 


; (K^(z) 

' 0 


* 

y 1 .h. (z))h, (z)dz  = 0 
. - , 11  k 

1 = 1 


wliich  is  an  inhomogeneous  set  of  linear  equations 


V 

k 


N 

T* 

/. 

i = l 


z,  . I . 
ki  1 


V = izii 


V 

k 


* 

K (z)li  (z)dz 

J ^ 

0 


'ki 


IL 

h. (z)h  (z)dz 

J 1 ^ 

0 


(9a) 


(9b) 
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with  solution 


I = [Z]  \ 


(10) 


a 

CO 

St 

fo 


Numerical  tests  on  the  above  algorithm  were  confined 
single  straight  wire  excited  across  a narrow  gap.  The 

llections  of  arbitrarily  oriented  skew  wires,  however, 
raightforward  since  the  integral  equation  analogous  to 
rm. 


to  cases  Involving 
extension  to 

is  reasonably 

(1  ) is  of  a similar 
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In  order  to  make  a comparison  between  the  MBCRK  alftorithm  and  other  MOM 

solutions  to  Equation  (1),  it  is  convenient  to  rewrite  (1)  in  operator  form. 

E (z)  = L(l)  (11a) 

z 

where  the  linear  operator  I.  is: 


L(  ) 


£ 

[ C(z,z' ) [ 
0 


]dz' 


(11b) 


If  the  inner  product  of  any  two  functions  along  the  wire  is  defined  as 

£ 


(f,g)  = 


f (z)g(z)dz 


(12) 


0 


then  the  error  of  equation  (5)  becomes 


= (r*,r) 


with 


r . - L(I^)  - E^  - 

The  generated  set  of  linear  equations  becomes 
* N ^ 

(h^.E^)  = [ h^)  k = 1,2, 


i=l 


or 


(L  I^L(f^))  = (L  (f^).r) 


. . N 


= 0 


(13a) 

(13b) 

(14a) 

(14b) 


Consider  now  an  arbitrary  MOM  solution  to  (1),  For  a given  set  of 
fj^(^),  the  solution  Is  completely  determined  by  the  choice  of  a set  of  M 
linearly  Independent  weighting  functions  Wj^(z).  The  generated  set  of 


linear  equations  becomes 


(15a) 


(ci.  , E ) 
k z 


(u),  , } I . L(f  . ) ) k = 1,2,  . . .M 

K , , A 1 

1=1 


or,  using  (13b) 


(u)j^,r)  = 0 


If  the  set  of  equations  (15b)  is  overdetermined , it  has  no  solution, 
but  a unique  pseudo-solution  can  be  obtained  by  defining  a discrete  residual 


error  vector  R and  a positive  definite  error  function  e,j,(R'  as 


r^  = 


r,  E R 
k — 


(16a) 


D V * 

®T  ° ^ '’k''k 

1 ,,=  1 k k 


The  minimization  of  e.^,  over  the  set  of  K is  analogous  to  the  minimization 


of  e^  in  Section  1.  Define 

I . = X.  + jy . 

J J J 


''k  = K’^z^ 


’^k  = ; ~jn  ^''k'.^,  =^ki^l^ 
1 1 


<1 


The  condition  for  a minimum  becomes 


V 1 * * 

^ 7^^ — ;T72 

J k=l 


= -2  Re[  I * ^--r/2  = 0 

k=l 


^®T  .V  1 * * 

° 7"* 

k=l 


Therefore 


k=l  ((^,1^,^) 


M 1 

7 * 1/2  *^k^kj  ° ° 


j = 1,  2,  ...  N 


(17b) 


In  matrix  form  this  can  be  written  as 


[Z]“[W]  = (Z]”[W]"^Z]I 


(18a) 


which  has  the  solution 


I = [[Z]“[W]"^Z]]'^[Z]”(W]"^V 


(18b) 


Here  [Z]  is  the  conjugate  transpose  of  [Z]  and  [W]  is  the  MxM  diagonal 
matrix  whose  elements  are  given  by 

In  this  way  any  moments  scheme  can  be  developed  from  an  error  minimiza- 
tion process.  Furthermore,  any  two  schemes  which  minimize  error  functions 
with  the  same  minimum  point  are  equivalent. 


If  the  set  of  weipliting  functions  is  taken  to  be  a contiguous  set  of 
unit  amplitude  pulses  of  equal  width  covering  the  wire,  the  inner  products 
of  ( 16a)  can  be  approx  limited  by 


Zj^+Az/2 


k * 1/2 

K’'“k^ 


r (z )dz 


r(Zk> 


z,  -Az/2 
k 


Az 

(Az) 


rz2 


(19) 


where  z is  the  coordinate  of  the  center  of  the  kth  pulse  and  Az  is  the 
k — 

11  . 

pulse  width.  The  error  function  e^  becomes 
M M 


I)  V * ^ ^ t \ r \ ^ 

T = Vk  = 

k=l  k=l 


(20) 


As  the  set  of  pulses  become  infinitely  dense,  the  finite  sum  passes  to  the 
defining  integral  for  e 


1 ■ 

lim  e = 
Az-0 


T 

‘k 


r (z)r(z)dz  = e_ 


(21) 


0 


Thus  the  MBCRE  algorithm  is  equivalent  to  a highly  overdetermined  point 
matching  formulation.  Conversely,  the  criteria  of  minimum  BCRE  can  be 
roughly  enforced  by  a variety  of  overdetermined  point  matching  schemes  in 
a way  that  is  numerically  cost  competitive  with  singly  determined  formulations. 

In  the  following  numerical  tests,  the  MBCRE  algorithm  is  compared 
with  two  other  algorithms  (pulse  weighting  and  plecewlse-slnusoidal  or 
c;alorkin  weighting)  all  using  piecewise-sinusoidal  expansion  functions. 

The  purpose  of  the  tests  is  to  determine  whether  the  MBCRE  algorithm  gives 
better  results  for  very  sparse  or  minimal  expansion  function  coverings  of 
thin  wires.  If  the  MBCRE  algorithm  were  significantly  better  under  these 
circumstances,  an  overdetermined  formulation  of  the  type  described  above 
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would  be  less  costly  to  apply  to  larqe  problems.  Unfortunately,  the 
MBCKK  algorithm  proved  to  be  generally  inferior  to  the  flalerkin  formulation, 
and  often  worse  than  the  pulse  formulation. 

F.xplanation  of  the  Tables 

The  following  set  of  test  problems  all  involve  a single  straight 

thin  wire  excited  across  a narrow  gap  by  a unit  voltage  source.  For  each 

test  problem,  a single  set  of  expansion  functions  f^  (z)  was  chosen,  along 

with  three  different  sets  of  weighting  functions  and  three  different 

functional  models  of  the  excitation  field  due  to  the  voltage  source.  The 

three  weighting  function  sets  (pulse,  Galerkin,  MBCRE)  were  defined  in 

a manner  shown  in  Figure  (lb).  The  MBCRE  weighting  functions  are 

(from  14b)  the  set  of  conjugate  fields  of  each  of  the  f.  (z) . The  three 

1 

gap  field  models  were  a unit  impulse  located  at  the  feed  point,  a one  volt 
pulse  centered  at  the  feed  point  with  a width  of  the  "gap  width"  of  the 
table,  and  a one  volt  isosceles  triangle  centered  at  the  feed  point  with 
a width  of  twice  the  gap  width  as  shown  in  Figure  (Ic).  For  each  test 
problem,  sets  of  amplitudes  were  calculated  for  all  possible 
combinations  of  weighting  functions  and  field  models,  making  a total  of 
nine  amplitude  sets  for  each  set  of  expansion  functions.  The  expansion 
functions  are  listed  with  their  three  defining  p>oints  z^,  and  z^  of 
Figure  (la). 
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The  tabulated  error  figures  "Datum  Relative  Error"  and  "Minimum  Possible 


Error"  (CRE  ) are  def  ined  below.  In  order  to  have  a measure  of  error  for 
m 

each  calculation  of  the  current  amplitudasl^  it  is  necessary  to  calculate  an 
approximate  "true  solution"  for  each  test  problem.  This  was  done  by 

the  method  of  moments  using  a very  dense  covering  of  piecewise  sinusoidal 
expansion  functions  with  Calerkin  weighting  and  using  a one  volt  impulse  to 
model  the  field  in  the  gap.  The  "true  solution"  was  verified  whenever  pos- 
sible by  comparison  with  measured  data  (Ref.  1)  or  with  results  obtained  by 
other  calculations  (Ref.  2).  Once  I,j,(z)  is  found,  it  is  possible  to  define 

for  each  set  of  I . a current  error  function  e(z)  and  a total  RMS  error  E by 
1 c ^ 


N 

e(z)  = I (z)  - I I ,f . (z) 

i=l  1 1 


(22a) 


E = (e  (z),e(z)) 
c 


(22b) 


The  set  of  currents  I . which  minimize  E are  derived  in  a manner  com- 
1 c 

pletely  analagous  with  the  minimization  of  e^  and  are  given  by  the  solution 
of  the  set  of  linear  equations 


(e(z),f.(z))  = 0 


i = 1,2 N 


The  error  function  can  now  be  broken  up  into  two  components 

N 

I 

i=l 


,(z)  = I„(z)  - I l°f,(z) 
i=l 

N N 

^(z)  = e(z)  - e^(z)  = I 


(23) 


(24a) 


(24b) 


Because  of  the  orthogonality  condition  (23): 


(e^(z),e2(z))  = (e^(z) ,62(2))  = 0 


(25) 


and  it  follows  that 
* 


(e  (z),e(z))  = (e^(z)  + 62 (2) ,e^ (z) + 62 (2) ) 


= (6j^(z)  ,e^^(z))  + (ej^(z)  ,62  (2))  + (e^(z),e2(z)) 


+ (02(2) ,62(2)) 


(26) 


|e^(z)|^  + |e2(z)|^ 


The  terms  "Datum  Relative  Error  (DRE)"  and  "Minimum  Possible  Error 


(CRE  )"  are  defined  as 
n 


CRE 


hj,  (z)  I 

m I I_(z) I 


(27a) 


DRE  = — 


62(2) I 


i-1 


i;fi(z) 


The  total  error  was  partitioned  in  this  way  because  CRE  is  an  error 

in 

which  is  due  entirely  to  the  choice  of  the  set  of  f^(z)  and  can  not  be  re- 
moved by  any  solution  algorithm.  The  DRE  is  the  error  component  which  is 
entirely  the  fault  of  the  solution  algorithm  and  is  therefore  a better 
measure  of  accuracy  than  the  total  RMS  error  E^. 
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Section  III 


The  computer  program  "MOP"  is  designed  to  be  run  on  the  Honeywell  GECOS 
time-sharing  system  in  the  Fortran  CARDIN  subsystem.  Input  data  is  read  from 
the  end  of  the  program  and  output  is  written  onto  a sequential  file  (with 
PRMFL  designation  "02"),  as  well  as  on  the  on-line  printer. 

The  following  is  a list  of  input  variables  and  their  units.  The  wire 
parameters  are 

WL  = wire  length  (meters) 

B - wave  number /2tt  = 1/A  (meters). 

RAD  = wire  radius  (meters) 

The  solution  for  the  I.j,(z)  is  specified  by: 

NPS  = number  of  expansion  functions  (always  equally  sized  and  equally 

spaced) 

KPS  = the  number  of  the  expansion  function  on  which  the  unit  impulse 
excitation  is  centered. 

The  test  problem  is  specified  by; 

NC  = Number  of  expansion  functions 

NE  = number  of  nodes  needed  to  define  the  three  expansion  function  end 
points  and  the  gap  end  points 
TE(K)  (K  = 1,NE)  = the  z coordinate  of  node  k. 

TI(K,.I)  (K  = 1,NC,  -J  = 1,3)=  Zj  (see  Eq,  3)  for  the  kth  expansion 
function 

KL  = node  number  for  one  end  of  the  gap  (lowest  z ). 

KH  = node  number  for  the  other  end  of  the  gap  (highest  z ) 

There  is  also  an  input  variable  KDEI.  which  is  the  number  of  the  expan- 


lA 


I 


1 


sion  function  amplitude  to  which  the  calculated  input  impedance  is  referred, 

i . e . 

Z.  = V.  /I  . = 1 ./KKDEI,) 

in  in  in 


The  calculation  of  Z.  will  only  be  meaningful  if  the  expansion  function 

in 

'kDEL  *^he  gap. 

All  variables  are  read  in  at  the  end  of  the  program  in  the  order  and 


with  the  formats  given  below: 
Wl.,  B,  RAD 
NFS,  KPS,  KDEL 
NC,  NE 

TE(K)  K = 1,NE 


(3F10.6) 

(314) 

(214) 

(F10.6) 

(each  node  on  a separate  I ine 1 


TI(K,J)  (K  = 1,NC),(J  = 1,3) 


(314) 

(The  three  defining  points  for  eacli 
expansion  function  on  one  line) 


KL,KH 


(214) 


Example : 

Assume  as  a test  problem  a center  excited  half  wave  dipole  of  radius 
O.OlA.  Assume  also  that  it  is  desired  to  test  the  three  MOM  formulations 
for  a set  of  three  equally  spaced  expansion  functions  of  equal  size  covering 
the  wire.  If  we  assume  an  excitation  gap  width  of  0.05A  and  a set  of  9 ex- 
pansion functions  for  l.j,(z),  the  data  cards  would  look  like  this. 
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<1 


The  program  output  is  written  onto  the  printer  and  specified  permanent 


file  in  the  following  order. 

(1)  The  input  impedance  calculated  from  I^(z) 

(2)  The  current  amplitudes  of  I.j,(z) 

(3)  The  test  problem  specifications: 

a)  wire  length 

b)  wave  number 

c)  wire  radius 

d)  the  three  defining  points  z z^,  z^  for  each  of  the  specified 
expansion  functions 

e)  the  end  points  of  the  gap. 

(A)  CRE 

m 

(5)  The  Input  Impedance  calculated  from  the  1° 

(6)  The  amplitudes  1° 

For  each  combination  of  excitation  field  model  and  weighting  function 
the  program  prints 

(7)  The  input  impedance  (referred  to 

(8)  The  DRE 

(9)  The 

Above  items  7,  8 and  9 will  appear  two  numbers  (LI  and  L2)  which 
both  vary  from  one  to  three  and  which  designate  the  type  of  excitation  field 
and  weighting  functions  l.e. 


LI  - 

1 

for  pulse  weighting 

LI  - 

2 

for  Galerkin  weighting 

LI  - 

3 

for  MBCRE  weighting 

L2  = 1 


for  impulse  excitation 


L2  = 2 for  pulse  excitation 

L2  = 3 for  triangle  excitation 

All  printed  currents  are  in  amperes  in  the  format 

K,  REAL  (1(K)),  AIMA(;(I(K)),  CMAC  (I(K)),  PHASE  (I(K)) 


References : 

1)  King,  Mack  and  Sandler,  Arrays  of  Cylindrical  Dipoles. 

Cambridge  University  Press,  1968. 

2)  Harrington  and  Mautz,  "Wires  with  Arbitrary  Excitation,"  IEEE  Trans. 
on  Antennas  and  Propagation,  July  1967. 


TEST  PROBLEM  ^ 


Wire  Length  0.50 


Radius  0.015 


Feed  Point  0.25 


Gap  Width  0.05 


Minimum  Possible 

Error  (CRE  ) 

0.049 

m *" 

Expansion  Functions 

No. 

"l 

*2 

*3 

1. 

0.000 

0.083 

0.166 

2. 

0.083 

0.166 

0.250 

3. 

0.166 

0.250 

0.333 

4. 

0.250 

0.333 

0.417 

5. 

0.333 

0.417 

0.500 

6. 

7. 

8. 

9. 

10. 

11. 

12. 
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TEST  PROBLEM  1 


Datum  Relative  Error 


Pulse 

Piecewise 

Sinusoid 

(Galerkin) 

MBCRE 

Impulse 

0.17 

0.045 

0.10 

Pulse 

0.17 

0.048 

0.13 

Triangle 

0.17 

0.049 

1 

0.15 

TEST  PROBLEM  2 


Wire  Length  p so Radius  0.015 


Feed  Point  0-  25  Gap  Width  0-125 


Minimum  Possible  Error  (CRE  ) 


Expansion  Functions 


No. 

"l 

^2 

^3 

1. 

0 , 000 

0.125 

0.250 

2. 

0.  125 

0.250 

0.375 

3. 

0.250 

0.375 

0.500 

A. 

5. 

6. 

7. 

8. 

9. 

10. 

11. 

12. 

22 


TEST  PROBLEM  ^ 


Datura  Relative  Error 


E^(^) 

Pulse 

Piecewise 

Sinusoid 

(Galerkin) 

MBCRE 

Impulse 

0.29 

0.078 

0.095 

Pulse 

0.29 

0.070 

0.33 

Triangle 

0.24 

0.091 

0.35 

TEST  PROBLEM  3 


Wire  Length  0»50 


Radius  0 . 015 


Feed  Point  0.25 


Gap  Width  0-05 


Minintua  Possible  Error  (CRE  ) 


Expansion  Functions 


Mo. 

*1 

*2 

*3 

1. 

0.000 

0.200 

0.250 

2. 

0.200 

0.250 

0.300 

3. 

0.250 

0.300 

0.500 

4. 

5. 

6. 

7. 

8. 

9, 

10. 

11. 

12. 


24 


TEST  PROBLEM  4 


Wire  Length  0. 50 


Radius  0.015 


Feed  Point  0.250 


Gap  Width  0.C5 


Minimum  Possible  Error  (CRE  ) 0.103 

m — — — 


Expansion  Functions 


Ho. 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9, 

10. 

11. 


0.000  0.250  0.500 

0.200  0.250  0.300 


12 
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TEST  PROBLEM 


Datum  Relative  Error 


Pulse 

i 

Piecewise 

Sinusoid 

(Galerkln) 

MBCRE 

Impulse 

0.78 

0.13 

0.42 

Pulse 

0.78 

0.14 

1 

1 

0.46 

Triangle 

0.77 

0.15 

0.47 

TEST  PROBLEM  5 


Wire  Length  ^ • OQO 


Radius  0.015 


Feed  Point 


0.5C0 


Gap  Width  0-100 


Minimum  Possible  Error  (CRE  ) O' 

m " 


Expansion  Functions 


No. 

^2 

*3 

1. 

0.000 

0.250 

0.400 

2. 

0.250 

0.400 

0.500 

3. 

0.  AOO 

0.500 

0.600 

4. 

0.  500 

0.600 

0.750 

s. 

0.600 

0.750 

1.000 

28 


. t 

I 


1 

TEST  PROBLEM  5 


Datum  Relative  Error 


i 

Pulse 

Piecewise 

Sinusoid 

(Galerkin) 

MBCRE 

Impulse 

0.13 

0.058 

0.32 

Pulse 

0,13 

0.125 

0.25 

Triangle 

0.22 

0.161 

0.28 

TEST  PROBLEM  6 


Wire  Length  1*0 


Radius  0.  C15 


Feed  Point  0.50 


Gap  Width  0-10 


Minimum  Possible  Error  (CRE  ) 

m 


0.090 


No. 

Expansion  Functions 

Zl  Z2 

*3 

1. 

0.000 

0.100 

0.200 

2. 

0.100 

0.200 

0.300 

3. 

0.200 

0.300 

0.400 

4. 

0.300 

0.400 

0.500 

5. 

0.400 

0.500 

0.600 

6. 

0.500 

0.600 

0.700 

7. 

0.600 

0.700 

0.800 

8. 

0.700 

0.800 

0.900 

0.800 

0.900 

1.000 

10. 

U. 


12 
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TEST  PROBLEM  ^ 


Datum  Relative  Error 


\w^(z) 

Pulse 

Piecewise 

Sinusoid 

(Galerkin) 

MBCRE 

Impulse 

0.  072 

0.020 

0.  336 

Pulse 

0.072 

0.097 

0.179 

Triangle 

0.165 

0.134 

0.212 

3 


I 


TEST  PROBLEM  ^ 


Wire  Length  1.0 


Radius  0.013 


Feed  Point  0.50 


Gap  Width  0.10 


Mlninum  Possible  Error 


(CRE^) 


0.141 


Expansion  Functions 

No. 

*1 

*2 

*3 

1. 

0.000 

0.166 

0.333 

2. 

0.166 

0.333 

0.500 

3. 

0.333 

0.500 

0.666 

4. 

0.500 

0.666 

0.833 

5. 

0.666 

0.833 

1.000 

6. 

7. 

8. 

9. 

10. 

11. 

12. 


<1 


TEST  PROBLEM 


7 


Datum  Relative  Error 


E (z) 
s 

Pulse 

Piecewise 

Sinusoid 

(Galerkin) 

MBCRE 

Impulse 

0.092 

0.045 

0.59 

Pulse 

0.092 

0.070 

0.18 

Triangle 

0.109 

0.105 

0.22 

33 


TEST  PROBLEM  8 


Wire  Length  !•  0 


Radius  0.015 


Feed  Point  0-50 


Gap  Width  0-10 


Minimum  Possible  Error  (CRE  ) 

m 


0.189 


Expansion  Functions 


^1 

^2 

*3 

0.000 

0.400 

0.450 

0.400 

0.450 

0.500 

0.450 

0.500 

0.550 

0.500 

0.550 

0.600 

0.550 

0.600 

1.000 

6. 

7. 

8. 

9, 

10. 

11. 

12. 


3A 


TEST  PROBLEM  8 


Datum  Relative  Error 


E^(z) 

Pulse 

Piecewise 

Sinusoid 

(Galerkin) 

MBCRE 

Impulse 

0.24 

0.29 

0.  27 

Pulse 

0.31 

0.29 

0.  30 

Triangle 

0.31 

0.33 

0.34 

35 


■4 


TEST  PROBLEM  9 


Wive  Length  1.0 


Radius  0.015 


Feed  Point  0.5 


Gap  Width  0.10 


Minimum  Possible  Error  (CRE  ) 

in 


0.19 


Expansion  Functions 


Mo. 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9, 

10. 

11. 


*1 

*2 

*3 

0.000 

0.250 

0.500 

0.400 

0.500 

0.600 

0.500 

0.750 

1.000 

12 
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TEST  PROBLEM  9 


Datum  Relative  Error 


Pulse 

Piecewise 

Sinusoid 

(Galerkln) 

1 

MBCRE 

Impulse 

0.42 

0.22 

0.36 

Pulse 

0.42 

0.31 

0.32 

Triangle 

0.54 

0.35 

1 

0.36 

TEST  PROBLEM 


Wire  Length  1-0  Radius  ..O-.OH 

Feed  Point  0-300  Gap  Width  _0^ 


Minimum  Possible  Error  (CRE  ) 

m 


0.20 


No. 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9, 

10. 

U. 


Expansion  Functions 


0.000 

0.400 


*2 

0.500 

0.500 


*3 

1.000 

0.600 


12 
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TEST  PROBLEM 


Datum  Relative  Error 


\wi(z) 

Pulse 

Piecewise 

Sinusoid 

(Galerkln) 

MBCRE 

Impulse 

0.22 

0.21 

0.36 

Pulse 

0.22 

0.31 

0.31 

Triangle 

1 

0.32 

0.35 

0.36 

TEST  PROBLEM 


Wire  Length  _ 

2.0 

Radius 

0.015 

Feed  Point 

l.O 

Gap  Width  ^ 

0.250 

Hlnlmuni  Possible  Error 

• (CRE  ) 

m 

0.13 

Expansion  Functions 

No. 

^1 

*2 

*3 

1. 

0.000 

0.250 

0.500 

2. 

1 

0.250 

0.500 

0.750 

1 

0.500 

0.750 

1.000 

■ 

. 4. 

0.750 

1.000 

1.250 

1 

i 5. 

1.000 

1.250 

1.500 

\ 

'1  6. 

1.250 

1.500 

1.750 

! 7. 

1.500 

1.750 

2.000 

8. 


11. 


12 


40 


TEST  PROBLEM 


11 


Datum  Relative  Error 


I 


I 


Pulse 

Piecewise 

Sinusoid 

(Galerkln) 

HBCRE 

Impulse 

0.13 

0.075 

0.75 

Pulse 

0.13 

0.185 

0.45 

Trlanp le 

0.33 

0.284 

0.51 
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